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I.  INTRODUCTION 


The  direct  coupling  approach  to  inertial  fusion,  in  which  an  imploding 
cylindrical  foil  liner  impacts  a  cylindrical  target  placed  inside  the  foil, 
is  unique  among  the  current  approaches  in  that  the  target  driving  energy  is 
already  in  the  form  of  inward  directed  kinetic  energy.  Since  a  goal  of 
fusion  target  design  is  to  convert  energy  of  a  given  initial  form  into 
directed  kinetic  energy  of  the  pusher  and  fuel,  the  fact  that  the  liner 
energy  is  already  in  inward  kinetic  energy  may  be  extremely  helpful.  No 
intermediate  energy  conversion  steps  are  necessary.  For  example,  in  fusion 
with  C02  lasers,  the  laser  energy  is  converted  into  the  energy  of  hot  elec¬ 
trons  at  the  critical  surface.  These  electrons  then  heat  the  ablator. 
Material  in  the  ablator  is  driven  off  at  high  velocities,  which  results  in 
inward  kinetic  energy  for  the  pusher  and  fuel.  In  this  process,  only  a 
small  fraction  of  the  laser  energy  is  ultimately  converted  into  kinetic 
energy  in  the  pusher  and  fuel.  Concepts  which  use  X-ray  radiation  to  drive 
targets  also  require  some  intermediate  energy  conversion  steps.  On  the 
other  hand,  the  cylindrical  symmetry  inherent  in  imploding  liner  systems, 
means  that  cylindrical  targets  must  probably  be  used  in  the  direct  coupling 
approach.  (Concepts  in  which  the  imploding  liner  is  coupled  to  a  spherical 
target  have  been  suggested,  but  their  feasibility  has  yet  to  be  demon¬ 
strated  (Refs.  1  through  3).)  This  presents  a  problem  since  it  is  more 
difficult  to  obtain  high  compression  in  a  cylindrical  target  than  in  a 
spherical  target. 

The  concept  of  imploding  a  cylindrical  deuterium-tritium  (DT)  target 
with  an  imploding  liner  has  received  considerable  attention  during  the  past 
decade.  However,  previous  work  has  focused  on  the  compression  and  heating 
of  low  density  DT  gas  by  a  relatively  slow  liner  (Refs.  4  through  12).  In 
general,  the  results  have  been  disappointing  in  that  very  high  input  ener¬ 
gies  were  required  to  achieve  only  modest  energy  gains.  Current  or  planned 
pulsed  powerexperiments  promise  higher  energies  with  shorter  implosion 
timescales.  Therefore,  it  Is  desirable  to  perform  a  study  with  input 
parameters  appropriate  to  this  new  generation  of  pulsed  power  machines. 
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The  SHIVA  experiment  at  the  Air  Force  Weapons  Laboratory  (AFWL)  is  of 
particular  interest  for  this  concept.  The  SHIVA  STAR  experiment,  which  is 
currently  being  made  operational,  has  a  6-MJ  fast  capacitor  bank  which  is 
used  with  inductive  storage--opening  switch  pulse  sharpening  to  drive  300 
to  500  ns,  30  to  50  cm/us,  400  kJ  to  1  MJ,  plasma  shells  having  axial 
length  on  the  order  of  2  cm.  Previous  experiments  at  AFWL  have  used  a  1.9 
MJ  capacitor  bank  to  implode  plasmas  with  kinetic  energy  200  to  250  KJ 
(Ref.  13). 

Thermonuclear  burn  in  spherical  targets  has  been  the  subject  of  intense 
study  during  the  past  10  years.  Consideration  of  the  Lawson  criterion  and 
numerous  numerical  studies  have  shown  that  bare  DT  spheres  must  be  com¬ 
pressed  so  that  /pdR  .>0.3  -  1.0  g/cm2  and  temperature  T  £  3-10  keV  in 
order  to  obtain  good  thermonuclear  burn.  In  addition.  Mason  and  Morse 
(Ref.  14)  have  shown  that  a  fixed  amount  of  fuel  gives  an  increased  yield 
for  a  given  p^R^  and  ignites  at  a  lower  p^  if  tamping  is  added,  where  pf 
and  are  the  fuel  density  and  radius. 

Obtaining  pfRf  >  0.3  g/cm2  is  likely  to  be  extremely  difficult  in 
cylindrical  systems,  as  shown  in  Section  II.  Hence,  it  is  important  to 
determine  if  significant  thermonuclear  burn  can  occur  in  tamped  cylindrical 
targets  which  have  pfRf  <  0.3  g/cm2.  Also,  an  axial  magnetic  field  can 
be  injected  into  a  cylindrical  system.  This  field  could  thermally  insulate 
the  DT  fuel  from  the  tamper  and  could  also  help  to  confine  the  alpha  par¬ 
ticles  emitted  by  the  0T  reactions.  Thus,  some  fraction  of  the  liberated 
energy  could  be  kept  in  the  fuel  instead  of  being  lost  to  the  tamper. 

In  Section  II,  the  difficulty  of  achieving  high  compression  in  cylin¬ 
drical  targets  is  illustrated.  In  Section  III,  estimates  of  the  magnetic 
field  strength  required  to  trap  the  alpha  particles  emitted  during  the  DT 
reaction  and  to  inhibit  the  electron  and  ion  thermal  conduction  are  made. 
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In  Section  IV,  calculations  of  thermonuclear  burn  in  tamped  targets  are 
presented.  The  effects  of  suppressing  the  thermal  conductivity  and 
trapping  the  alpha  particles  are  considered.  In  Section  V,  low  compression 
targets  are  considered  and  the  conclusion  is  reached  that  these  targets  are 
not  appropriate  for  fast  liner  systems  such  as  the  SHIVA  system.  In  Sec¬ 
tion  VI,  models  for  the  production  of  large  fluxes  of  neutrons  with  the 
SHIVA  STAR  experiment  are  considered. 

The  major  conclusion  from  this  study  is  that  cylindrical  targets  are 
not  suitable  targets  for  fast  liner  experiments  like  the  SHIVA  STAR.  The 
compression  requirements  for  the  fuel  are  too  severe,  given  the  symmetry  of 
fast  imploding  plasma  devices.  However,  large  quantities  of  neutrons  (1015 
to  3  x  1016)  may  be  produced  in  the  SHIVA  system.  This  could  be  extremely 
useful  for  diagnostic  purposes. 


II.  COMPRESSION  REQUIREMENTS 


The  purpose  of  this  section  is  to  illustrate  the  difficulty  in  obtain¬ 
ing  high  compression  in  cylindrical  targets  as  opposed  to  the  relative  ease 
in  spherical  targets.  As  mentioned  in  the  Introduction,  the  Lawson  cri¬ 
terion  and  numerous  numerical  studies  show  that  / pdr  =  0.3  to  1.0  g/cm2  for 
good  thermonuclear  burn  in  untamped  targets.  Clearly,  it  is  easier  to 
obtain  higher  /pdr  in  targets  with  larger  fuel  mass.  However,  the  energy 
requirements  are  such  that  only  a  fraction  of  a  milligram  of  DT  can  be 
brought  to  burn  temperature. 


For  a  constant  density  sphere  of  mass  1115,  density  ps,  and  radius 


R 


s» 


p 


sRs 


(i: 


For  a  corresponding  cylinder  of  length  L, 


M* 

PA  = 


(2) 


In  a  sphere  pR  a  1/R2,  but  in  a  cylinder  pR  a  1/R^.  Consider  an  example 
where  pR  =  1  g/cm2  and  m  =  0. 1  mg.  For  a  constant  density  sphere.  Equa¬ 
tion  1  gives  a  density  of  205  g/cm2,  which  is  a  factor  103  higher  density 
than  solid  0T.  For  a  corresponding  cylinder  with  L  =  1  cm.  Equation  2 
gives  a  density  of  3.14  x  104  g/cm3,  which  is  a  factor  of  more  than  105 
higher  density  than  solid  DT.  In  spherical  systems  it  is  experimentally 
difficult  to  achieve  compressions  of  103  over  solid  density;  compressions 
of  a  factor  10s  are  clearly  impossible  given  the  symmetry  of  current  avail 
able  in  Imploding  linear  systems. 


Thus,  systems  in  which  massive  tampers  hold  the  DT  fuel  together  are 
studied.  In  such  systems,  only  volume  ignition  of  the  fuel  is  possible;  it 
is  not  realistic  to  expect  hot  spot  Ignition  and  propagating  burn  as  will 
be  shown  in  the  next  section. 
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III.  MAGNETIC  FIELDS  TO  INHIBIT  THERMAL  CONDUCTION  AND  TRAP  ALPHA 

PARTICLES 


The  idea  that  magnetic  fields  can  be  used  to  inhibit  thermal  conduction 
and  trap  the  alpha  particles  emitted  during  the  DT  reactions  is  not  new. 
This  idea  has  been  examined  by  numerous  researchers  (Ref.  15)  over  the 
years.  This  section  presents  simple  estimates  of  the  magnetic  field 
strength  required  to  be  effective. 

Stronger  magnetic  fields  are  required  to  inhibit  the  ion  thermal  con¬ 
duction  than  to  inhibit  the  electron  thermal  conduction.  Also,  in  some 
cases  of  interest,  the  electron  and  ion  temperatures  will  not  be  equal.  It 
is  of  interest  to  examine  ion  and  electron  conduction  separately. 


The  magnetic  field  can  inhibit  thermal  conduction  in  a  direction  per¬ 
pendicular  to  the  magnetic  field  when  we  xe  >  1,  where  is  the 
electron  cyclotron  frequency  given  by 

ojgCys’1)  =  1.76  x  107  B(MG)  (3) 


and  x 0  is  the  electron  self-collision  time  given  by 


Te(M$) 


1.835  x 


A  Te372 
T  pfcnA 


(4) 


Te  is  the  electron  temperature  in  keV,  p  is  the  density  in  g/cm3,  A  is 
the  atomic  weight,  Z  is  the  atomic  number,  and  tnA  is  the  coulomb  logorithm 
(Refs.  16,  17).  Furthermore,  the  electron  thermal  conductivity  will  be 
reduced  over  the  classical  value  (without  magnetic  fields)  by  a  factor  of 
approximately  Fe  =  (u)eTe)2  if  wefe  >  1.  So 


F 

e 


B2(MG)T;J 

0.104 - £ 

p 


(5) 


There  is  a  strong  dependence  on  B,  T,  and  p.  Thus,  the  thermal  conducti¬ 
vity  can  more  easily  be  suppressed  in  a  low  density,  low  compression  target 


r 


than  in  a  high  density,  high  compression  target.  For  A/2  =2.5  appropriate 
for  OT  and  tnA  =  5,  Equation  5  becomes 


F  ■ 
e 


-2 

p 


(6) 


which  shows  that  the  electron  thermal  conduction  can  be  reduced  by  a  huge 
factor  in  imploding  liner  devices. 


The  ion  thermal  conductivity  is  inhibited  by  the  magnetic  field  when 
tu-j t -j  >  1,  where  u-j  is  the  ion  cyclotron  frequency  given  by 


^(ps-1)  =  9.6  x  103  (fy  B(MG) 
and  ri  is  the  ion  self  collision  time  given  by 


A3/2t3/2 

Mus)  =  1.113  x  1(T6  — r— ! - 

Z  ptn\ 


Defining  F-j  =  (oi-j t-j  )2  , 


F.  =  1.4  x  10 


-4  e2(MG>Tf  » 


"7  l6(in\)Z 


or 


( 8  f  (  t<  y 

i  -  vrocy  \m~wjj 


-2 


(7) 


(8) 


(9) 


where  Z  =  1,  A  =  2.5,  and  tnA  =  5.  Equation  9  shows  that  the  ion  thermal 
conduction  can  also  be  inhibited  in  imploding  liner  devices  if  the  density 
is  not  too  high.  However,  stronger  magnetic  fields  are  required  to  inhibit 
the  ion  conduction  than  to  inhibit  the  electron  conduction  for  equal  ion 
and  electron  temperatures. 


For  <ut  »  1,  the  thermal  conductivity  is  reduced  by  a  factor  of  approxi¬ 
mately  (<*>t)2  over  the  value  with  no  magnetic  field.  However,  Lindemuth  and 
Widner  (Ref.  18)  point  out  that  when  wt  *  1,  a  degree  of  magnetothermoi n- 
sulation  can  be  provided.  A  wete  of  unity  leads  to  a  reduction  factor 
of  nearly  4. 
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In  unmagnetized  fuel,  the  ratio  of  the  alpha  particle  range  \  to  the 
fuel  radius  Rf  is  given  by  approximately 


Thus,  in  a  high  compression,  untamped  target,  the  confinement  time  require¬ 
ment  (Lawson  criterion)  of  pRf  >  0.3  to  1.0  g/cm2  is  usually  more  severe 
than  the  alpha  trapping  requi rements. 

For  the  magnetic  field  to  aid  in  trapping  the  3.5  MeV  alpha  particles, 
the  gyroradius  of  the  alpha  particle  Ra  must  usually  be  smaller  than  the 
fuel  radius.  (One  could  imagine  situations  where  this  is  not  the  case. 

For  example,  the  fuel  region  could  be  surrounded  by  a  magnetic  field.  As 
long  as  the  alpha  particles  did  not  reach  the  wall  or  tamper,  but  only 
remained  in  the  vacuum  region  or  in  the  fuel,  the  alphas  would  eventually 
deposit  their  kinetic  energy  in  the  fuel  region.)  The  radius  of  gyration 
for  a  3.5  MeV  alpha  particle  is  given  by 

Ra  (cm)  =  fiMsy  (11) 

Thus,  it  will  be  extremely  difficult  to  trap  the  alpha  particles  in  a  high 
compression  (small  fuel  radius)  target  by  means  of  a  magnetic  field.  In  a 
low  compression  target  where  the  fuel  is  held  together  by  a  tamper  or  wall, 
the  magnetic  field  can  be  of  considerable  benefit. 


IV.  TAMPED  CYLINDRICAL  TARGETS 


It  was  shown  in  Section  II,  that  obtaining  compressions  in  cylindrical 
targets  so  that  / pR  a  0.3  to  1.0  g/cm2  will  be  extremely  difficult.  The 
degree  to  which  a  tamper  can  reduce  the  compression  requirements  is  inves¬ 
tigated  in  this  section  along  with  the  effects  of  reducing  the  thermal  con¬ 
ductivity  and  trapping  of  the  alpha  particles  as  might  be  expected  in  the 
presence  of  a  sufficiently  strong  magnetic  field. 

1.  DESCRIPTION  OF  CALCULATIONS 

The  computations  were  performed  with  the  one-dimensional  (1-D)  magneto¬ 
hydrodynamic  radiation  code  MACH1  (Ref.  19  and  20).  The  initial  conditions 
for  the  fuel  and  tamper  were  chosen  as  follows.  The  central  fuel  region 
was  a  cylinder  of  DT  at  initial  density  pf  and  initial  temperature  Tf. 

This  was  surrounded  by  a  nickel  shell  (the  tamper)  at  initial  density 
and  initial  temperature  Tt.  As  the  fuel  will  have  been  compressed  by  the 
tamper,  which  also  serves  as  a  pusher,  the  pressure  in  the  tamper  is 
required  to  be  equal  to  that  in  the  fuel.  For  all  the  computations 
reported  here,  the  fuel  mass  per  unit  length  was  chosen  to  be  0.1  mg/cm. 
This  fuel  mass  is  appropriate  for  a  SHIVA-1  ike  device  because  it  takes  58 
kJ  of  energy  to  heat  0.1  mg  of  DT  to  5  keV.  It  is  convenient  to  use  the 
variable  p^  as  an  independent  variable  since  the  burn  fraction  is 
expected  to  be  a  strong  function  of  p^.  IT  the  mass  per  unit  length  is 
fixed,  p^Rf  determines  p^  and  R^.  For  most  of  the  calculations  reported 
here  pt  is  assumed  to  equal  4  pf,  although  at  the  end  of  Section  III,  the 
effects  of  varying  pfc  are  discussed.  Choosing  a  value  of  T^,  then  gives  Tt 
if  pt  is  known.  Hence,  the  input  parameters  are  PfRf>  tamper  mass/length 
mt>  Tf,  and  pt«  Figure  1  is  a  schematic  of  the  initial  density  distribu¬ 
tion. 


Radiative  losses  were  not  included  for  the  calculations  reported  here. 
The  MACH1  code  runs  much  faster  if  radiation  is  not  included.  In  limited 


Tamped  Cylindrical  Targets 
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Figure  1.  Schematic  of  the  initial  density  distribution  used  in  calculations. 

tests,  it  was  found  that  radiative  losses  were  not  important  for  the 
parameter  range  of  most  interest.  In  any  case,  radiative  losses  will  only 
make  the  computed  thermonuclear  burn  smaller.  Classical  thermal  conducti¬ 
vity  was  either  included  or  suppressed  totally  as  one  might  expect  in  a 
strong  magnetic  field.  The  alpha  particles  emitted  during  the  DT  reaction 
were  assumed  to  either  escape  the  fuel  or  to  be  trapped  locally  as  might  be 
expected  in  the  presence  of  a  sufficiently  strong  magnetic  field. 

2.  RESULTS 

Figures  2  and  3  show  the  fraction  of  the  DT  fuel  burned  as  a  function 
of  tamper  mass  for  p^R^  =0.1  g/cm2  and  =  0.01  g/cm2,  respectively. 

For  all  points  plotted  in  Figures  2  and  3,  pt  =  4  pf. 
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Figure  3.  Fraction  of  DT  fuel  burned  versus  tamper 
mass  for  =  0.01  g/cm^. 


a.  pfRf  =0.1  g/cm2 --The  open  circles  in  Figure  2  are  results  for 
models  with  Tf  =  3  keV,  local  trapping  of  the  a  particles,  and  classical 
thermal  conductivity.  The  burn  fraction  increases  with  tamper  mass  from 
1.3  x  10-4  for  a  bare  DT  sphere  to  0.035  for  a  tamper  mass  of  1000  mf. 

This  increase  is  almost  solely  due  to  the  increased  confinement  time  as  the 
tamper  mass  is  increased.  Trapping  of  the  alpha  particles  is  only  impor¬ 
tant  in  the  m^  =  1000  mf  case.  In  models  with  lower  tamper  masses,  the 
heating  due  to  the  alpha  particles  is  not  enough  to  raise  the  temperature 
significantly.  A  reasonable  tamper  mass  for  a  SHIVA-like  device  is  about 
102  mf,  so  Tf  =  3  keV  and  p^R^  =0.1  g/cm2  are  not  sufficient  for  good 
burn. 

The  filled  circles  in  Figure  2  are  results  for  models  identical  to 
those  described  above  except  that  Tf  =  6  keV.  These  cases  are  much  more 
promising  since  tamper  masses  greater  than  100  mf  give  burn  fractions  in 
excess  of  10  percent.  Trapping  of  the  alpha  particles  is  critical,  how¬ 
ever.  The  filled  squares  are  results  when  the  alpha  particle  heating  is 
suppressed.  The  burn  fraction  is  only  2  percent  at  m^  =  100  mf. 

Trapping  of  the  alpha  particles  is  not  required  if  the  DT  fuel  can  be 
brought  to  a  sufficiently  high  temperature.  The  filled  triangles  in  Figure 
2  are  results  for  Tf  =  10  keV  and  no  heating  by  the  alpha  particles.  The 
open  triangles  are  results  when  the  thermal  conductivity  is  suppressed. 

Burn  fractions  approaching  10  percent  are  obtained  for  m^  =  100  mf. 

Higher  burn  fractions  would  be  obtained  for  higher  initial  temperatures. 

b.  pfRf  =  0.01  g/cm2--Fi gure  3  shows  the  burn  fraction  as  a  func¬ 
tion  of  tamper  mass  for  pfRf  =  0.01  g/cm2.  The  results  are  much  worse 
than  in  the  pfRf  =0.1  g/cm2  cases.  The  open  circles  are  results  for 
models  with  Tf  =  3  keV,  local  trapping  of  the  a  particles,  and  classical 
thermal  conductivity.  The  burn  fraction  ranges  from  1.2  x  10-5  for  a  bare 
cylinder  to  2.7  x  10‘4  for  a  tamper  with  1000  mf.  These  values  are  at 
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least  an  order  of  magnitude  worse  than  the  corresponding  models  with 
pfRf  =0.1  g/cm2. 


The  filled  circles  are  results  where  Tf  =  6  keV  and  include  classical 
thermal  conductivity.  The  results  are  an  order  of  magnitude  higher  than  in 
the  3  keV  models,  but  still  too  low  for  good  burn.  In  both  the  3  and  6  keV 
models,  it  makes  little  difference  whether  the  alpha  particles  are  trapped 
or  not.  There  is  not  enough  burn  to  increase  the  gas  temperature  signifi¬ 
cantly. 

Some  increase  in  burn  is  gained  by  suppressing  the  thermal  conductivi¬ 
ty.  This  is  illustrated  by  the  triangles  in  Figure  3,  which  are  results 
for  Tf  =  6  keV  with  suppressed  thermal  conductivity.  The  gain  is  not 
sufficient,  however,  to  give  good  burn  as  the  burn  fraction  is  only  3.8  x 
10“3  at  mf  =  100  mf. 

The  squares  in  Figure  3  are  results  for  extremely  optimistic  condi¬ 
tions.  In  these  models  Tf  =  10  keV,  the  thermal  conductivity  is  totally 
suppressed,  and  the  alpha  particles  are  assumed  to  be  trapped.  The  burn 
fraction  is  only  10-2  at  mf  =  100  mf. 

Hence  for  pfRf  <  0.01,  there  appears  to  be  no  reasonable  way  to 
achieve  good  burn  in  tamped  cylindrical  targets. 

c.  Thermal  Conducti vity--In  these  previous  results,  it  was  assumed 
that  pt  =  4  pf.  One  might  wonder  if  the  burn  fraction  could  be 
increased  by  using  more  dense  tampers  or  more  sophisticated  tampers.  The 
answer  is  no.  Figure  4  shows  the  burn  fraction  as  a  function  of  tamper 
density  for  pfRf  =  0.01  g/cm2,  Tf  =  6  keV,  and  m^  =  1000  mf.  The  burn 
fraction  actually  decreases  dramatically  as  the  tamper  density  increases. 
This  decrease  is  due  almost  entirely  to  thermal  conductivity.  As  the 
tamper  density  is  increased,  the  tamper  temperature  is  lowered  to  maintain 
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pressure  equilibrium.  The  denser  tampers  cause  excessive  energy  loss  from 
the  fuel  to  the  tamper  via  thermal  conduction.  In  the  case  where  pf  = 
pf,  the  tamper  is  hotter  than  the  fuel  and  energy  is  actually  transferred 
to  the  fuel  via  thermal  conduction.  Low  density  tampers  are  not  reasonable 
for  target  designs.  As  the  tamper  density  is  decreased,  the  temperature 
must  be  increased  to  maintain  pressure  equilibrium,  and  the  energy  require¬ 
ments  become  prohibitive.  For  example,  the  model  with  =  pf  gives  a 
burn  fraction  of  0.0261,  but  requires  76  MJ/cm  in  internal  energy  to 
achieve  the  initial  configuration.  Clearly,  such  models  are  not 
reasonable. 

When  thermal  conduction  was  totally  suppressed  in  models  with  pR  =  0.01 
g/cm2,  it  was  found  that  the  burn  fraction  was  almost  independent  of  tamper 
density  but  depended  strongly  on  tamper  mass.  Hence,  the  burn  fractions 
shown  in  Figure  3  where  the  thermal  conduction  is  totally  suppressed  are 
almost  independent  of  tamper  density.  Since  the  compressed  state  of  the 
tamper-pusher  and  fuel  will  be  in  approximate  pressure  equilibrium,  it  is 
likely  that  the  use  of  complex  tamper  will  not  significantly  increase  the 
burn  fractions  above  those  calculated  above  as  long  as  the  energy  require¬ 
ments  are  reasonable.  While  this  result  has  not  been  proven  in  any 
rigorous  sense,  as  it  is  based  on  computer  models  which  sample  a  limited 
range  of  parameter  space,  we  feel  that  it  is  valid  for  the  parameter  range 
of  most  interest. 

3.  CONCLUSIONS  ANO  DISCUSSION 

The  goal  in  this  section  has  been  to  determine  the  conditions  in  the  0T 
fuel  and  tamper  In  order  to  achieve  significant  thermonuclear  burn  in 
cylindrical  systems.  If  the  fuel  can  be  compressed  to  pfRf  =  1.0 
g/cm2,  then  good  burn  can  be  obtained  if  Tf  >  3  keV.  For  pfRf  =0.1 
g/cm2,  reasonably  good  burn  can  be  obtained  if  the  alpha  particles  are 
trapped,  a  tamper  of  100  mf  or  larger  is  used,  and  Tf  >  6  keV. 
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Reasonable  burn  can  also  be  obtained  if  Tf  >  10  keV,  even  if  the  alpha 
particles  are  not  trapped.  For  pfRf  =  0.01  g/cm2,  significant  burn 
could  not  be  achieved  for  reasonable  parameters,  even  with  the  most  opti¬ 
mistic  assumptions.  These  were  that  the  alpha  particles  were  trapped, 
thermal  conduction  was  totally  suppressed,  and  Tf  =  10  keV. 

To  this  point  the  energy  requirements  for  the  models  presented  have  not 
been  discussed  in  any  detail.  For  massive  tampers  the  energy  requirements 
are  critically  dependent  on  tamper  density.  The  energy  required  to  raise 
the  DT  fuel  is  easily  calculated  from  the  temperature 

E(MJ)  =  0.115  mf(mg)  Tf(kev) 

For  example,  it  takes  115  kJ  to  raise  0.1  mg  of  DT  to  10  keV.  The  tamper 
is  likely  to  be  in  pressure  equilibrium  with  the  fuel.  So,  the  denser  the 
tamper  can  be  made,  the  less  energy  is  required.  For  example,  the  energy 
required  for  the  models  discussed  in  connection  with  Figure  4  ranged  from 
0.42  MJ/cm  to  76  MJ/cm  as  the  tamper  density  decreased  from  120  pf  to 

1  Pf. 


From  the  target  design  standpoint,  the  tamper  must  be  made  as  dense  as 
possible  in  order  to  reduce  the  energy  requirements.  At  high  tamper  densi 
ties,  thermal  conduction  must  be  suppressed  so  that  the  energy  in  the  fuel 
is  not  lost  to  the  tamper. 

Reluctantly,  it  must  be  concluded  that,  unless  targets  can  be  com¬ 
pressed  so  that  pfRf  -  0.1  g/cm2,  then  a  good  burn  (>  10  percent) 
cannot  be  achieved  in  tamped  cylindrical  systems  with  reasonable 
parameters.  For  a  0.1  mg  target,  pfRf  =0.1  g/cm2  requires  a  fuel  den¬ 
sity  of  314  g/cm3,  a  compression  of  greater  than  1000  over  solid  density. 
Given  the  symmetry  In  current  imploding  liner  devices,  this  compression  is 
extremely  unlikely  to  be  achieved  in  the  near  future. 
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V.  LOW  COMPRESSION  TARGETS 


It  was  shown  in  the  previous  section  that  it  is  difficult  to  achieve 
good  thermonuclear  burn  in  tamped  targets  unless  high  fuel  compression  can 
be  achieved.  In  this  section,  a  slightly  different  concept  is  evaluated  in 
which  the  OT  fuel  is  held  together  by  the  ram  pressure  of  an  imploding 
plasma  liner.  The  model  presented  is  very  simple  and  crude,  but  it  is 
sufficient  to  show  the  difficulties  of  this  approach  for  fast  liner  systems 
such  as  SHIVA.  More  detailed  models  are  described  in  References  4  and  5 
and  7  through  10. 

Suppose  that  some  DT  gas  has  been  brought  to  thermonuclear  burn  temper¬ 
atures  (-10  keV)  in  some  unspecified  manner.  It  is  necessary  to  estimate 
the  energy  required  for  a  fast  plasma  liner  to  hold  this  DT  gas  together 
long  enough  for  breakeven  (so  that  the  thermonuclear  burn  energy  is  equal 
to  the  energy  required  to  bring  the  DT  gas  to  10  keV).  It  is  assumed  that 
there  is  sufficient  magnetic  field  at  the  liner-fuel  interface  so  that 
there  are  no  thermal  conduction  losses. 

The  Lawson  criterion  for  breakeven  is 

Nt  *  1014  s/cm3  (12) 

where  N  is  the  number  density  of  hydrogen  nucleii  and  t  is  the  confinement 
time.  The  pressure  P  in  the  DT  gas  is  given  by 

P  =  2NkT  =  2  x  1014  (13) 

where  it  has  been  assumed  that  the  ion  and  electron  temperatures  T  are 
equal.  Boltzmann's  constant  is  k.  In  order  to  confine  the  plasma,  the  ram 
pressure  in  the  liner  must  be  equal  to  the  gas  pressure  in  the  DT  fuel. 
Then, 
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(14) 


d  2  2K 

P  =  pv  =  — 7 
irl< 


where  p,  v,  K  are  the  liner  density,  velocity,  and  kinetic  energy  per  unit 
length.  R  is  the  radius  of  the  cylindrical  liner.  It  has  been  assumed 
that  the  fuel  radius  is  much  smaller  than  the  liner  radius.  Combining 
Equations  13  and  14, 

K  =  2  x  1014  irR2kT/t  (15) 


Making  the  approximation  t  =  R/v,  Equation  15  can  be  written  as 

K’5Hj/c"fc)  N  W  U6) 

For  a  fast  liner  system  such  as  SHIVA,  v  -  40  cm/ps,  R  -  1  cm.  The 
energy  required  to  confine  the  plasma  for  breakeven  is  clearly  prohibitive 
for  fast  liner  systems. 


From  Equation  16,  K  is  proportional  to  vR.  Hence,  if  a  DT  plasma  is  to 
be  combined  with  a  cylindrical  liner,  a  slow  liner  which  remains  thin 
(perhaps  even  remains  solid)  is  required  so  that  the  energy  requirements 
are  not  enormous.  This  concept  has  been  suggested  by  a  group  at  Los  Alamos 
National  Laboratory  (Ref.  10).  Another  approach  would  be  to  confine  the 
plasma  with  rigid  walls,  supply  enough  magnetic  field  for  sufficient 
magnetothermoinsulation,  and  shock  heat  the  plasma  by  some  other  means  as 
suggested  by  Kmetyk  and  Gross  (Refs.  21  and  22). 
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VI.  NEUTRON  PRODUCTION 


1.  INTRODUCTION  AND  MOTIVATION 

The  previous  sections  have  shown  that  good  thermonuclear  burn  in  cylin¬ 
drical  targets  is  extremely  difficult  with  fast  imploding  plasma  liners 
such  as  the  SHIVA  STAR  system.  This  section  shows  that  a  simple  system  can 
produce  large  amounts  of  neutrons  (1015  -  3  x  1016)  whicn  may  be  useful  for 
diagnostic  purposes. 

The  concept  put  forward  here  is  as  follows.  Using  gas  puff  techniques 
place  a  small  mass  (about  1/100  the  mass  of  the  liner)  in  a  column  along 
the  axis  of  symmetry.  When  the  fast  plasma  liner  impacts  this  low  density 
gas,  it  is  shock  and  compressional ly  heated  to  thermonuclear  burn  tempera¬ 
ture.  There  are  probably  many  different  concepts  which  one  could  be  used 
to  produce  large  amounts  of  neutrons  in  an  imploding  liner  device.  Some 
advantages  of  the  proposed  concept  are  as  follows.  The  concept  and  calcu¬ 
lations  based  on  the  concept  are  very  simple  and,  hence,  are  more  reliable 
than  a  complex  design.  Gas  puff  techniques  have  been  used  in  the  SHIVA 
program  before  (Refs.  23  and  24).  Hence,  the  technology  already  exists  at 
AFWL  to  perform  the  experiment.  Large  neutron  production  does  not  require 
good  symmetry  in  the  imploding  plasma  liner.  Even  if  hot  spots  are  formed 
in  the  target,  significant  neutrons  will  be  emitted  from  these  hot  spots. 
Lastly,  if  the  emitted  neutrons  can  be  imaged,  they  can  be  a  powerful  diag¬ 
nostic  in  determining  the  quality  and  symmetry  of  the  plasma  implosion. 

The  most  likely  reason  for  the  possible  failure  of  this  concept  would  be 
substantial  mixing  between  the  liner  material  and  DT  gas  before  thermonu¬ 
clear  burn  conditions  are  achieved  in  the  gas. 

Previously,  liners  which  have  some  deuterium  in  them  have  been  used  to 
produce  thermonuclear  neutrons  at  AFWL.  This  is  not  a  good  method  if  the 
quality  of  the  plasma  implosion  is  good.  One-dimensional  calculations  show 
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that  the  neutron  production  in  such  systems  is  very  low.  In  particular, 
one  calculation  was  performed  where  a  0.5  cm  thick  plasma  of  pure  DT  with 
an  energy  of  350  kJ/cm  and  a  velocity  of  40  cm/ys  was  allowed  to  pinch  on 
axis.  The  neutron  production  in  this  model  was  only  7  x  1012  n/cm.  In 
this  calculation,  the  maximum  ion  temperature  was  less  than  1  keV  because 
the  internal  energy  was  rapidly  radiated  away.  If  some  amount  of  higher  Z 
material  had  been  added  to  the  DT,  the  ion  temperature,  and  hence,  the  neu¬ 
tron  yield  would  have  been  much  lower  since  the  plasma  would  have  been  able 
to  radiate  much  more  effectively.  However,  in  two-dimensional  calculations 
or  in  experiments,  the  situation  may  be  somewhat  different  because  insta¬ 
bilities  can  cause  the  formation  of  low  density  regions  (which  may  also 
have  higher  velocity)  which  can  be  shock  heated  and  compressed  by  the  main 
portion  of  the  imploding  plasma.  Such  a  neutron  production  method  relies 
on  a  poor  quality  implosion.  It  is  more  desirable  to  inject  the  low  den¬ 
sity,  pure  DT  gas  into  the  experiment. 

2.  DESCRIPTION  OF  CALCULATIONS  AND  SELECTION  OF  PARAMETERS 

The  computations  were  performed  with  the  1-D  magnetohydrodynamic  radia¬ 
tion  code  MACH1  (Refs.  19  and  20). 

The  following  simple  model  was  chosen  for  the  computations.  An  outer 
plasma  liner  of  uniform  density,  velocity,  and  temperature  impacts  a  cylin¬ 
drical  DT  target  of  uniform  density  and  temperature  at  time  t  =  0.  The 
computation  is  continued  at  least  until  the  neutron  production  rate  falls 
to  a  negligible  level.  Parameters  appropriate  to  the  SHIVA  STAR  experiment 
were  chosen  for  the  outer  plasma  sheath.  A  typical  calculation  has  initial 
kinetic  energy  E  =  350  kJ/cm,  velocity  v  =  40  cm/ys,  thickness  aRl  =  0.5 
cm,  and  is  assumed  to  be  aluminum.  For  the  DT  fuel,  a  typical  mass  was 
0.05  mg,  and  a  typical  fuel  radius  Rf  was  0.5  cm.  All  of  the  above 
parameters  were  varied  within  reasonable  limits.  In  all  models  the  initial 
temperature  of  the  plasma  sheath  was  5  eV,  and  the  initial  temperature  of 
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the  DT  fuel  was  1  eV.  Input  parameters  for  the  33  models  reported  on  here 
are  listed  in  Table  1. 

Effects  due  to  radiation  were  found  to  be  important,  so  radiation  was 
included  in  all  models  except  two. 

Thermal  conductivity  was  also  found  to  be  extremely  important.  The 
electron  thermal  conduction  is  especially  important  at  the  fuel-liner 
interface.  As  shown  in  Section  II,  the  electron  thermal  conductivity  will 
be  inhibited  due  to  strong  magnetic  fields.  It  is  difficult  to  reliably 
compute  the  magnetic  field  strength  near  this  interface.  The  theta  com¬ 
ponent  of  the  magnetic  field  Bg  depends  on  the  current  flowing  interior 
to  the  interface  which  will  be  difficult  to  compute  during  the  pinch  phase 
of  the  implosion.  Even  if  the  Bg  field  is  not  strong  enough  to 
sufficiently  impede  the  electron  conductivity,  a  small  Bz  field  ( <10  kG) 
injected  prior  to  the  start  of  the  liner  implosion  can  provide  the  required 
magnetothermoinsulation  at  the  interface.  In  the  calculations  presented 
here,  the  electron  thermal  conductivity  is  either  totally  suppressed  or  it 
is  allowed  to  flow  as  if  there  were  no  magnetic  field.  In  Table  1,  the 
column  labeled  t.c.  shows  if  electron  thermal  conductivity  is  allowed. 

The  ion  thermal  conductivity  is  not  so  important  at  the  fuel-liner 
interface.  The  conductivity  coefficient  is  usually  small  in  the  liner 
material  because  of  the  1/Z4  dependence  (see  Eq.  8).  Ion  conduction  is 
important  at  the  center  of  the  fuel  region.  Without  ion  conduction  a  cen¬ 
tral  hot  spot  forms  in  the  innermost  zone.  However,  it  will  be  difficult 
to  inhibit  the  ion  conduction  at  the  center  of  the  fuel  region.  By 
symmetry  the  Bg  field  vanishes  at  the  center,  and,  hopefully,  a  strong 
Bz  field  will  not  have  diffused  to  the  center  of  the  fuel  or  it  will 
impede  the  fuel  compression.  For  these  reasons,  the  ion  conduction  is  not 
inhibited  in  the  calculations  presented  in  this  report. 
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TABLE  1.  NEUTRON  PRODUCTION 


Fuel 


Model 

No. 

E 

(kJ/cm) 

V 

cm/us 

mass 

mg 

Rf 

cm 

aR  [_ 
cm 

t.c. 

RAD 

Material 

Yield 

(1015/cm 

1 

350 

40 

0.1 

1/2 

1/2 

N* 

Y* 

At 

2.1 

2 

350 

40 

0.05 

1/2 

1/2 

N 

Y 

At 

2.6 

3 

350 

40 

0.02 

1/2 

1/2 

N 

Y 

At 

2.1 

4 

350 

40 

0.01 

1/2 

1/2 

N 

Y 

At 

1.3 

5 

350 

40 

0.1 

1/2 

1/2 

Y 

Y 

At 

0.036 

6 

350 

40 

0.05 

1/2 

1/2 

Y 

Y 

At 

0.089 

7 

350 

40 

0.02 

1/2 

1/2 

Y 

Y 

At 

0.60 

8 

350 

40 

0.01 

1/2 

1/2 

Y 

Y 

At 

1.1 

9 

350 

40 

0.005 

1/2 

1/2 

Y 

Y 

At 

0.87 

10 

350 

40 

0.05 

1/2 

1/2 

N 

Y 

Ni 

4.0 

11 

350 

40 

0.05 

1/2 

1/2 

N 

Y 

CH 

1.6 

12 

350 

40 

0.05 

1/2 

1/2 

N 

Y 

Be 

2.2 

13 

350 

40 

0.05 

1/2 

1/2 

Y 

Y 

Ni 

0.051 

14 

350 

40 

0.05 

1/2 

1/2 

Y 

Y 

CH 

0.13 

15 

350 

40 

0.05 

1/2 

1/2 

Y 

Y 

Be 

0.22 

16 

350 

40 

0.05 

1/2 

1/2 

N 

N 

At 

4.6 

17 

350 

40 

0.05 

1/2 

1/2 

Y 

N 

At 

0.73 

18 

200 

40 

0.02857 

1/2 

1/2 

N 

Y 

At 

1.0 

19 

500 

40 

0.07143 

1/2 

1/2 

N 

Y 

At 

4.9 

20 

200 

40 

0.02857 

1/2 

1/2 

Y 

Y 

At 

0.10 

21 

500 

40 

0.07143 

1/2 

1/2 

Y 

Y 

At 

0.091 

22 

350 

40 

0.05 

1 

1/2 

N 

Y 

At 

8.1 

23 

350 

40 

0.05 

2 

1/2 

N 

Y 

At 

9.8 

24 

350 

40 

0.05 

1 

1 

N 

Y 

At 

1.7 

25 

350 

40 

0. 05 

1 

1/2 

Y 

Y 

At 

0.98 

26 

350 

40 

0.05 

2 

1/2 

Y 

Y 

At 

4.6 

27 

350 

40 

0.05 

1 

1 

Y 

Y 

At 

0.18 

28 

350 

40 

0.1 

2 

1/2 

N 

Y 

At 

15.3 

29 

350 

40 

0.02 

2 

1/2 

N 

Y 

At 

4.3 

30 

350 

40 

0.1 

2 

1/2 

Y 

Y 

At 

2.0 

31 

350 

40 

0.02 

2 

1/2 

Y 

Y 

At 

4.0 

32 

350 

20 

0.05 

1/2 

1/2 

N 

Y 

At 

1.2 

33 

350 

20 

0.05 

1/2 

1/2 

Y 

Y 

At  6. 

5  x  10-4 

*N  *  No,  Y  =  Yes 
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A  final  point  is  that  some  care  has  been  exercised  in  insure  that  the 
MACH1  code  treats  both  the  electron  and  ion  conductivity  correctly  in 
regions  where  there  are  large  temperature  gradients  (Ref.  25). 

3.  RESULTS 

Although  the  major  result  of  the  calculations,  the  neutron  yield,  is 
listed  in  Table  1,  the  results  from  two  models  are  considered  first  in  some 
detail.  Model  6  is  a  case  with  the  typical  parameters  described  above  and 
where  the  electron  thermal  conductivity  is  not  suppressed.  Figures  5 
through  11  show  the  electron  and  ion  temperatures,  the  energy  flux  due  to 
electron  thermal  conduction  and  to  radiation,  the  velocity,  and  the  density 
plotted  versus  radius  at  various  times  for  Model  6.  Figure  5  shows  these 
quantities  at  time  t  =  0  as  discussed  above  (the  initial  conditions). 

Figure  6  shows  that  the  fuel -liner  interface  has  moved  to  a  radius  of 
0.36  cm.  A  strong  shock  is  being  driven  into  the  fuel  region.  The  energy 
dissipated  in  the  shock  is  transferred  to  the  ions,  and  the  ion  temperature 
is  higher  than  the  electron  temperature  just  behind  the  shock.  Electron 
thermal  conduction  is  carrying  energy  ahead  of  the  shock  so  the  electron 
temperture  is  higher  than  the  ion  temperature  near  a  radius  of  0.30  cm. 

Note  that  a  single  strong  shock  brings  the  ion  temperature  in  the  fuel 
region  to  only  0.6  keV — not  enough  to  achieve  good  burn.  The  density  in 
the  fuel  region  is  low.  The  density  in  the  liner  shows  a  sharp  spike  near 
the  interface.  The  energy  flux  plot  shows  how  energy  is  loss  from  the  fuel 
region.  Electron  thermal  conduction  transports  energy  from  the  fuel  to  the 
edge  of  the  liner  where  it  is  radiated  away.  The  fuel  region  itself  is  not 
dense  enough  to  radiate  effectively,  but  the  liner  is.  It  might  be 
expected  that  either  shutting  off  the  thermal  conduction  or  reducing  the 
ability  of  the  liner  to  radiate  would  decrease  the  loss  of  energy  in  the 
fuel  and  thereby  increase  the  neutron  yield. 
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Figure  6.  Temperature,  energy  flux  velocity,  and  density  at  time  t  = 
for  Model  6.  Tj  is  ion  temperature.  Te  is  electron  temper 
ature.  TC  is  electron  thermal  flux.  R  is  radiation  flux. 


:ure  (keV) 


rature  (keV) 


rature  (keV) 


Figure  7  shows  that  the  interface  is  at  0.22  cm  at  t  =  8  ns.  The  first 
strong  shock  has  almost  reached  the  center. 

At  t  =  12  ns,  the  interface  is  at  0.08  cm.  The  ion  temperature  in  the 
fuel  is  2  keV.  The  fuel  density  has  increased  by  a  factor  30  over  the  ini¬ 
tial  value.  The  thermal  conduction  is  very  large  reaching  a  peak  value  of 
1.7  x  1012  W. 

At  t  ■  14  ns,  the  interface  is  at  0.03  cm.  The  ion  temperature  is  in 
excess  of  3  keV.  Energy  losses  due  to  electron  thermal  conduction  are 
enormous.  This  is  the  time  of  maximum  thermonuclear  burn. 

At  t  =  16  ns,  the  ion  temperature  has  fallen  to  1  keV.  Energy  is 
rapidly  being  transported  from  the  center.  Since  the  temperature  cannot  be 
maintained  in  the  fuel  region,  there  is  not  enough  pressure  to  keep  the 
liner  from  compressing  the  fuel  to  high  density. 

At  t  =  18  ns,  the  fuel  temperature  has  dropped  to  less  than  0.3  keV, 
and  the  thermonuclear  burn  is  essentially  over,  while  the  plasma  liner  is 
still  collapsing  on  axis. 

Figure  12  shows  the  ion  temperature,  the  density,  the  pressure,  and  the 
velocity  of  zones  1  (the  central  zone),  10  (half  way  through  the  fuel 
region  by  mass),  20  (outer  edge  of  fuel),  21  (inner  edge  of  plasma  liner), 
30  (a  point  in  the  interior  of  the  liner)  as  a  function  of  time.  The  peak 
in  the  Ion  temperature  occurs  before  the  peak  in  the  density.  The  thermo¬ 
nuclear  yield  rate  and  the  energies  as  a  function  of  time  are  also  shown  in 
Figure  12.  The  yield  rate  has  a  sharp  rise  and  an  even  sharper  fall.  Most 
of  the  original  kinetic  energy  is  ultimately  converted  to  radiation.  After 
a  time  of  10  ns,  there  is  approximate  pressure  equilibrium  in  the  fuel  and 
the  Inner  region  of  the  liner.  In  Section  IV,  a  pressure  equilibrium  was 
assumed  between  the  fuel  and  tamper.  The  pressure  curves  Illustrate  the 
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Figure  12.  Ion  temperature,  density,  pressure,  velocity  yield  rate,  and 
energy  histories  for  Model  6.  In  the  upper  four  plots,  the 
quantities  are  for  zones  I  (central  zone),  10  (half  way  through 
the  0T  fuel  region  by  mass),  20  (outer  edge  of  fuel  region),  21 
(inner  edge  of  liner),  and  30  (in  the  liner,  initially  at 
radius  5.3  cm).  In  the  energy  plot,  KE  is  the  kinetic  energy, 

I  the  internal  energy,  and  RAO  is  the  energy  radiated  away. 


reasonableness  of  this  assumption.  There  is  some  evidence  for  multiple 
shocks  in  the  velocity  plot  for  zone  10. 

Model  2  is  identical  to  Model  6  except  that  the  electron  thermal  con¬ 
duction  is  totally  suppressed.  Figure  5  is  a  plot  of  the  initial  condi¬ 
tions  for  Model  2  as  well  as  for  Model  6.  Figures  13  through  18  for  Model 
2  should  be  compared  to  Figures  6-10  for  Model  6. 

Figure  13  is  very  similar  to  Figure  6  at  time  4  ns.  The  temperatures 
in  the  shocked  region  are  higher  in  Model  2  since  the  electron  thermal  con¬ 
duction  is  suppressed.  The  radiation  flux  is  slightly  less  in  Model  2. 

Figures  14  and  15  for  Model  2  should  be  compared  to  Figures  6  and  8  for 
Model  6.  The  curves  are  again  similar.  The  temperatures  are  higher  and 
the  radiation  flux  is  lower  in  Model  2. 

However,  when  Figure  16  is  compared  to  Figure  10  there  is  a  large 
difference.  The  fuel  region  in  model  2  remains  hot— there  is  no  electron 
thermal  conductivity  to  cool  it.  Since  it  is  not  cooled,  the  fuel  pressure 
remains  high,  and  the  fuel  density  stays  comparatively  low.  The  minimum 
interface  radius  is  0.035  cm,  which  is  a  fuel  compression  of  a  factor  200. 
This  is  the  time  of  maximum  thermonuclear  burn. 

Figures  17  and  18  show  that  the  fuel  region  remains  hot  even  as  it 
starts  to  expand.  The  ion  temperature  is  still  4  keV  at  time  28  ns. 

Figure  19  shows  histories  for  Model  2  which  should  be  compared  to  the 
histories  for  Model  6  shown  in  Figure  12.  The  fuel  ion  temperature  is 
higher  and  the  curve  is  broader  for  Model  2.  The  fuel  density  reaches  a 
lower  peak  in  Model  2.  The  velocity  curve  for  zone  10  clearly  shows  tht 
multiple  shocks  have  been  driven  through  the  fuel.  The  maximum  internal 
energy  is  higher  in  Model  2. 
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Ion  temperature,  density,  pressure,  velocity,  yield  rate,  and 
energy  histories  for  Model  2.  In  the  upper  four  plots,  the 
quantities  are  for  zones  1  (central  zone),  10  (half  way  through 
the  DT  fuel  region  by  mass),  20  (outer  edge  of  fuel  region),  21 
(Inner  edge  of  liner),  and  30  (In  the  liner.  Initially  at 
radius  5.3  cm).  In  the  energy  plot,  KE  Is  the  kinetic  energy, 

I  the  internal  energy,  and  RAO  Is  the  energy  radiated  away. 


The  effects  of  variation  of  the  input  parameters  are  now  considered. 

The  last  column  of  Table  1  gives  the  neutron  production  for  the  33  models 
reported  on  here. 

In  Models  1  through  9,  the  fuel  mass  was  varied.  Models  1  through  4, 
in  which  the  electron  thermal  conduction  is  suppressed,  show  little  varia¬ 
tion  in  neutron  yield.  The  maximum  production  of  2.6  x  1015  n/cm  is 
reached  for  a  fuel  mass  of  0.05  mg  (Model  2).  Models  5  through  9,  in  which 
the  conduction  is  not  suppressed,  show  a  large  variation  in  neutron  yield. 
The  maximum  production  of  1.1  x  1015  n/cm  is  reached  for  a  fuel  mass  of 
0.01  mg.  Lower  fuel  masses  are  required  if  the  conduction  cannot  be 
suppressed  by  magnetic  fields. 

Models  2,  6,  and  10  through  15  show  the  effects  of  different  liner 
materials.  When  the  conduction  is  suppressed,  the  variation  in  neutron 
production  is  small  although  high  Z  liners  give  somewhat  better  yields.  If 
the  conduction  cannot  be  suppressed  then,  lower  Z  liners  give  better 
yields.  This  is  easily  understood,  since  the  region  in  the  liner  near  the 
fuel -liner  interface  cannot  radiate  as  effectively  in  low  Z  materials.  In 
Models  16  and  17  the  radiation  is  turned  off.  These  models  are  of  no 
practical  interest  since  there  is  no  way  to  turn  the  radiation  off  in  an 
experiment. 

Models  2,  6,  and  18  through  21  give  the  neutron  production  when  the 
liner  kinetic  energy  is  changed.  When  the  conduction  is  suppressed,  the 
yield  is  strongly  dependent  on  the  input  energy,  varying  from  1.0  x  1015  to 
4.9  x  1015  n/cm  as  the  kinetic  energy  is  varied  from  200  to  500  kJ/cm. 

There  is  very  little  variation  in  production  when  the  electron  conduction 
Is  not  suppressed. 

Models  2,  6,  and  22  through  27  give  the  neutron  production  when  the 
thickness  of  the  liner  and  the  thickness  of  the  fuel  region  are  varied. 
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Clearly,  thicker  fuel  regions  give  much  higher  neutron  yields.  A  thicker 
liner  results  in  poorer  performance.  The  fuel  mass  in  Models  28  through  31 
is  varied  as  well  as  the  thickness.  Model  28  had  the  highest  neutron  pro¬ 
duction  (1.5  x  1016  n/cm)  of  any  model. 

Models  32  and  33  show  what  happens  when  the  liner  velocity  is  lowered 
to  20  cm/ps.  When  conduction  is  suppressed,  the  yield  is  a  factor  2.2 
lower  than  in  Model  2.  When  conduction  is  not  suppressed,  the  yield  is  a 
factor  140  lower  than  in  Model  6. 

In  conclusion,  one  would  like  to  chose  thick  (>2  cm  radius)  DT  fuel 
regions  and  one  would  like  to  keep  the  plasma  liner  as  thin  as  possible. 
Fuel  masses  of  about  0.1  mg  are  best  if  one  believes  that  the  electron  con¬ 
duction  can  be  suppressed  (as  is  most  probably  the  case).  The  liner  is 
about  2  cm  long  in  the  SHIVA  STAR  system.  So  the  results  from  Model  28 
would  give  a  neutron  yield  of  3  x  1016  n  for  this  sytem.  When  electron 
thermal  conduction  is  supressed,  there  is  suprisingly  little  variation  in 
the  neutron  yield  as  the  input  parameters  are  varied.  Indeed,  most  models 
in  Table  1  with  supressed  conduction,  could  be  roughly  characterized  as 
giving  a  neutron  yield  of  a  few  x  1015/cm.  This  lack  of  variation  in  yield 
gives  one  considerable  confidence  that  an  experiment  on  the  SHIVA  STAR 
system  would  be  successful. 
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VII.  FINAL  COMMENTS 


It  has  been  shown  that  it  will  be  extremely  difficult  to  achieve  good 
thermonuclear  burn  (>10  percent)  in  cylindrical  targets  driven  by  fast 
plasma  liners  with  realistic  parameters.  Tamped  cylindrical  target  systems 
require  high  compression  (p^R^  *  0.1  g/cm2)  to  achieve  good  burn.  Given 
the  symmetry  in  the  current  generation  of  imploding  liner  devices,  this 
compression  is  extremely  unlikely  to  be  achieved  in  the  near  future.  For 
low  compression,  wall  confined  targets  with  magnetothermoinsulation,  other 
methods  may  be  better  choices. 

A  large  neutron  flux  (1015  -  3  x  1016)  can  be  produced  in  a  simple  fast 
liner-target  system  appropriate  for  the  SHIVA  STAR  experiment.  The  neutron 
yield  was  not  very  sensitive  to  the  input  parameters  as  long  as  the  elec¬ 
tron  thermal  conductivity  was  suppressed  at  the  liner  fuel  interface.  The 
neutron  production  was  highest  when  the  fuel  radius  was  large  (>2  cm). 
Hussey  (Ref.  26)  has  suggested  that  a  low  density  gas  fill  may  help  to 
stabilize  the  plasma  liner  during  the  run  in  phase  of  the  implosion.  If 
the  fill  gas  is  DT,  then  the  neutrons  produced  could  be  a  powerful  diagnos¬ 
tic  in  determining  the  quality  of  the  plasma  implosion.  In  the  future,  it 
would  be  desirable  to  have  one-dimensional  and  two-dimensional  magneto- 
hydrodynamic  calculations  of  this  concept. 
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